Lessons from prospective longitudinal follow-up of a French APECED cohort.

BACKGROUND: APECED syndrome is a rare disease caused by biallelic mutations of the AIRE gene, usually presenting with the triad "hypoparathyroidism-adrenal failure-chronic mucocutaneous candidiasis (CMC)" and non-endocrine manifestations. The aim of this study was to determine the molecular profile of the AIRE gene, the prevalence of rare manifestations and to characterize immunological disturbances in a French cohort. PATIENTS AND METHODS: A national, multicenter prospective observational study to collect genetic, clinical, biological and immunological data (NCT03751683). RESULTS: 25 patients (23 families) were enrolled. Eleven distinct AIRE variants were identified, two of which were not previously reported: an intronic variant, c.653-70G > A, and a c.1066del (p.Arg356GlyfsX22) variant (exon 9). The most common was the Finnish variant c.769C > T (16 alleles), followed by the variant c.967_979del13 (15 alleles), which seemed associated with a less severe phenotype. 17/25 patients were homozygote. The median number of clinical manifestations was seven; 19/25 patients presented with the hypoparathyroidism-adrenal failure-CMC triad, 8/13 showed pulmonary involvement, 20/25 had ectodermal dystrophy, 8/25 had malabsorption, and 6/23 had asplenia. Fifteen out of 19 patients had NK cell lymphopenia with an increase in CD4+ and CD8+ T lymphocytes and an age-dependent alteration of B lymphocyte homeostasis compared with matched controls (p < 0.001), related to the severity of the disease. All tested sera (n = 18) were positive for anti-interferon-α, 15/18 for anti-interleukin-22 antibodies, and 13/18 for anti-interleukin-17F antibodies, without clear phenotypic correlation other than with CMC. CONCLUSION: This first prospective cohort showed a high AIRE genotype variability, with two new gene variants. The prevalence of potentially life-threatening non-endocrine manifestations, was higher with systematic screening. These manifestations could, along with age-dependent B-cell lymphopenia, contribute to disease severity. Systematic screening for all the manifestations of the syndrome would allow earlier diagnosis, supporting vaccination, and targeted therapeutic approaches.

Heterozygous loss-of-function variants in DOCK4 cause neurodevelopmental delay and microcephaly.

Neurons form the basic anatomical and functional structure of the nervous system, and defects in neuronal differentiation or formation of neurites are associated with various psychiatric and neurodevelopmental disorders. Dynamic changes in the cytoskeleton are essential for this process, which is, inter alia, controlled by the dedicator of cytokinesis 4 (DOCK4) through the activation of RAC1. Here, we clinically describe 7 individuals (6 males and one female) with variants in DOCK4 and overlapping phenotype of mild to severe global developmental delay. Additional symptoms include coordination or gait abnormalities, microcephaly, nonspecific brain malformations, hypotonia and seizures. Four individuals carry missense variants (three of them detected de novo) and three individuals carry null variants (two of them maternally inherited). Molecular modeling of the heterozygous missense variants suggests that the majority of them affect the globular structure of DOCK4. In vitro functional expression studies in transfected Neuro-2A cells showed that all missense variants impaired neurite outgrowth. Furthermore, Dock4 knockout Neuro-2A cells also exhibited defects in promoting neurite outgrowth. Our results, including clinical, molecular and functional data, suggest that loss-of-function variants in DOCK4 probable cause a variable spectrum of a novel neurodevelopmental disorder with microcephaly.

Assessment of parental mosaicism rates in neurodevelopmental disorders caused by apparent de novo pathogenic variants using deep sequencing.

While de novo variants (DNV) are overall at low risk of recurrence in subsequent pregnancies, a subset is at high risk due to parental mosaicism. Accurately identifying cases of parental mosaicism is therefore important for genetic counseling in clinical care. Some studies have investigated the rate of parental mosaics, but most were either limited by the sensitivity of the techniques (i.e. exome or genome sequencing), or focused on specific types of disease such as epileptic syndromes. This study aimed to determine the proportion of parental mosaicism among the DNV causing neurodevelopmental disorders (NDDs) in a series not enriched in epilepsy syndromes. We collected 189 patients with NDD-associated DNV. We applied a smMIP enrichment method and sequenced parental blood DNA samples to an average depth of 7000x. Power simulation indicated that mosaicism with an allelic fraction of 0.5% would have been detected for 87% of positions with 90% power. We observed seven parental mosaic variants (3.7% of families), of which four (2.1% of families) had an allelic fraction of less than 1%. In total, our study identifies a relatively low proportion of parental mosaicism in NDD-associated DNVs and raises the question of a biological mechanism behind the higher rates of parental mosaicism detected in other studies, particularly those focusing on epileptic syndromes.

Penetrance, variable expressivity and monogenic neurodevelopmental disorders.

PURPOSE: Incomplete penetrance is observed for most monogenic diseases. However, for neurodevelopmental disorders, the interpretation of single and multi-nucleotide variants (SNV/MNVs) is usually based on the paradigm of complete penetrance. METHOD: From 2020 to 2022, we proposed a collaboration study with the French molecular diagnosis for intellectual disability network. The aim was to recruit families for whom the index case, diagnosed with a neurodevelopmental disorder, was carrying a pathogenic or likely pathogenic variant for an OMIM morbid gene and inherited from an asymptomatic parent. Grandparents were analyzed when available for segregation study. RESULTS: We identified 12 patients affected by a monogenic neurodevelopmental disorder caused by likely pathogenic or pathogenic variant (SNV/MNV) inherited from an asymptomatic parent. These genes were usually associated with de novo variants. The patients carried different variants (1 splice-site variant, 4 nonsense and 7 frameshift) in 11 genes: CAMTA1, MBD5, KMT2C, KMT2E, ZMIZ1, MN1, NDUFB11, CUL3, MED13, ARID2 and RERE. Grandparents have been tested in 6 families, and each time the variant was confirmed de novo in the healthy carrier parent. CONCLUSION: Incomplete penetrance for SNV and MNV in neurodevelopmental disorders might be more frequent than previously thought. This point is crucial to consider for interpretation of variants, family investigation, genetic counseling, and prenatal diagnosis. Molecular mechanisms underlying this incomplete penetrance still need to be identified.

High diagnostic potential of short and long read genome sequencing with transcriptome analysis in exome-negative developmental disorders.

Exome sequencing (ES) has become the method of choice for diagnosing rare diseases, while the availability of short-read genome sequencing (SR-GS) in a medical setting is increasing. In addition, new sequencing technologies, such as long-read genome sequencing (LR-GS) and transcriptome sequencing, are being increasingly used. However, the contribution of these techniques compared to widely used ES is not well established, particularly in regards to the analysis of non-coding regions. In a pilot study of five probands affected by an undiagnosed neurodevelopmental disorder, we performed trio-based short-read GS and long-read GS as well as case-only peripheral blood transcriptome sequencing. We identified three new genetic diagnoses, none of which affected the coding regions. More specifically, LR-GS identified a balanced inversion in NSD1, highlighting a rare mechanism of Sotos syndrome. SR-GS identified a homozygous deep intronic variant of KLHL7 resulting in a neoexon inclusion, and a de novo mosaic intronic 22-bp deletion in KMT2D, leading to the diagnosis of Perching and Kabuki syndromes, respectively. All three variants had a significant effect on the transcriptome, which showed decreased gene expression, mono-allelic expression and splicing defects, respectively, further validating the effect of these variants. Overall, in undiagnosed patients, the combination of short and long read GS allowed the detection of cryptic variations not or barely detectable by ES, making it a highly sensitive method at the cost of more complex bioinformatics approaches. Transcriptome sequencing is a valuable complement for the functional validation of variations, particularly in the non-coding genome.

Congenital hydrocephalus: new Mendelian mutations and evidence for oligogenic inheritance.

BACKGROUND: Congenital hydrocephalus is characterized by ventriculomegaly, defined as a dilatation of cerebral ventricles, and thought to be due to impaired cerebrospinal fluid (CSF) homeostasis. Primary congenital hydrocephalus is a subset of cases with prenatal onset and absence of another primary cause, e.g., brain hemorrhage. Published series report a Mendelian cause in only a minority of cases. In this study, we analyzed exome data of PCH patients in search of novel causal genes and addressed the possibility of an underlying oligogenic mode of inheritance for PCH. MATERIALS AND METHODS: We sequenced the exome in 28 unrelated probands with PCH, 12 of whom from families with at least two affected siblings and 9 of whom consanguineous, thereby increasing the contribution of genetic causes. Patient exome data were first analyzed for rare (MAF < 0.005) transmitted or de novo variants. Population stratification of unrelated PCH patients and controls was determined by principle component analysis, and outliers identified using Mahalanobis distance 5% as cutoff. Patient and control exome data for genes biologically related to cilia (SYScilia database) were analyzed by mutation burden test. RESULTS: In 18% of probands, we identify a causal (pathogenic or likely pathogenic) variant of a known hydrocephalus gene, including genes for postnatal, syndromic hydrocephalus, not previously reported in isolated PCH. In a further 11%, we identify mutations in novel candidate genes. Through mutation burden tests, we demonstrate a significant burden of genetic variants in genes coding for proteins of the primary cilium in PCH patients compared to controls. CONCLUSION: Our study confirms the low contribution of Mendelian mutations in PCH and reports PCH as a phenotypic presentation of some known genes known for syndromic, postnatal hydrocephalus. Furthermore, this study identifies novel Mendelian candidate genes, and provides evidence for oligogenic inheritance implicating primary cilia in PCH.

Retrotransposon insertion as a novel mutational cause of spinal muscular atrophy.

Spinal muscular atrophy (SMA) is an autosomal recessive neuromuscular disorder resulting from biallelic alterations of the SMN1 gene: deletion, gene conversion or, in rare cases, intragenic variants. The disease severity is mainly influenced by the copy number of SMN2, a nearly identical gene, which produces only low amounts of full-length (FL) mRNA. Here we describe the first example of retrotransposon insertion as a pathogenic SMN1 mutational event. The 50-year-old patient is clinically affected by SMA type III with a diagnostic odyssey spanning nearly 30 years. Despite a mild disease course, he carries a single SMN2 copy. Using Exome Sequencing and Sanger sequencing, we characterized a SINE-VNTR-Alu (SVA) type F retrotransposon inserted in SMN1 intron 7. Using RT-PCR and RNASeq experiments on lymphoblastoid cell lines, we documented the dramatic decrease of FL transcript production in the patient compared to subjects with the same SMN1 and SMN2 copy number, thus validating the pathogenicity of this SVA insertion. We described the mutant FL-SMN1-SVA transcript characterized by exon extension and showed that it is subject to degradation by nonsense-mediated mRNA decay. The stability of the SMN-SVA protein may explain the mild course of the disease. This observation exemplifies the role of retrotransposons in human genetic disorders.

Variants in CLDN5 cause a syndrome characterized by seizures, microcephaly and brain calcifications.

The blood-brain barrier ensures CNS homeostasis and protection from injury. Claudin-5 (CLDN5), an important component of tight junctions, is critical for the integrity of the blood-brain barrier. We have identified de novo heterozygous missense variants in CLDN5 in 15 unrelated patients who presented with a shared constellation of features including developmental delay, seizures (primarily infantile onset focal epilepsy), microcephaly and a recognizable pattern of pontine atrophy and brain calcifications. All variants clustered in one subregion/domain of the CLDN5 gene and the recurrent variants demonstrate genotype-phenotype correlations. We modelled both patient variants and loss of function alleles in the zebrafish to show that the variants analogous to those in patients probably result in a novel aberrant function in CLDN5. In total, human patient and zebrafish data provide parallel evidence that pathogenic sequence variants in CLDN5 cause a novel neurodevelopmental disorder involving disruption of the blood-brain barrier and impaired neuronal function.

Contribution of DNA methylation profiling to the reclassification of a variant of uncertain significance in the KDM5C gene.

KDM5C encodes a demethylase of the histone H3 lysine 4 residue, involved in chromatin regulation and gene expression. Hemizygous KDM5C pathogenic variants cause X-linked intellectual disability of Claes-Jensen type. Because of its mode of inheritance and the low specificity of the clinical phenotype, interpretation of variants can be difficult, hence the need for functional studies and biomarkers specific to this disorder. We present the case of a male patient with intellectual disability, behavioral abnormalities and subtle dysmorphic features, in which genetic investigation identified a hemizygous novel missense KDM5C variant of uncertain significance (VUS), inherited from his asymptomatic mother and present in his paucisymptomatic sister. We assessed the global genomic DNA methylation status from a whole blood sample of the proband. Global DNA methylation profiling specifically identified the recently discovered epi-signature of Claes-Jensen syndrome. This result served as a biomarker which independently highlighted KDM5C as the cause of the disorder in this patient. Because of the X-linked mode of inheritance, variant reclassification had a high impact on genetic counseling in this family. This example highlights the value of global methylome profiling in situations of variants of uncertain significance in genes with a known specific epi-signature.

Deep intronic NIPBL de novo mutations and differential diagnoses revealed by whole genome and RNA sequencing in Cornelia de Lange syndrome patients.

Cornelia de Lange syndrome (CdLS; MIM# 122470) is a rare developmental disorder. Pathogenic variants in 5 genes explain approximately 50% cases, leaving the other 50% unsolved. We performed whole genome sequencing (WGS) ± RNA sequencing (RNA-seq) in 5 unsolved trios fulfilling the following criteria: (i) clinical diagnosis of classic CdLS, (ii) negative gene panel sequencing from blood and saliva-isolated DNA, (iii) unaffected parents' DNA samples available and (iv) proband's blood-isolated RNA available. A pathogenic de novo mutation (DNM) was observed in a CdLS differential diagnosis gene in 3/5 patients, namely POU3F3, SPEN, and TAF1. In the other two, we identified two distinct deep intronic DNM in NIPBL predicted to create a novel splice site. RT-PCRs and RNA-Seq showed aberrant transcripts leading to the creation of a novel frameshift exon. Our findings suggest the relevance of WGS in unsolved suspected CdLS cases and that deep intronic variants may account for a proportion of them.

Toward clinical and molecular dissection of frontonasal dysplasia with facial skin polyps: From Pai syndrome to differential diagnosis through a series of 27 patients.

Unique or multiple congenital facial skin polyps are features of several rare syndromes, from the most well-known Pai syndrome (PS), to the less recognized oculoauriculofrontonasal syndrome (OAFNS), encephalocraniocutaneous lipomatosis (ECCL), or Sakoda complex (SC). We set up a research project aiming to identify the molecular bases of PS. We reviewed 27 individuals presenting with a syndromic frontonasal polyp and initially referred for PS. Based on strict clinical classification criteria, we could confirm only nine (33%) typical and two (7%) atypical PS individuals. The remaining ones were either OAFNS (11/27-41%) or presenting with an overlapping syndrome (5/27-19%). Because of the phenotypic overlap between these entities, OAFNS, ECCL, and SC can be either considered as differential diagnosis of PS or part of the same spectrum. Exome and/or genome sequencing from blood DNA in 12 patients and from affected tissue in one patient failed to identify any replication in candidate genes. Taken together, our data suggest that conventional approaches routinely utilized for the identification of molecular etiologies responsible for Mendelian disorders are inconclusive. Future studies on affected tissues and multiomics studies will thus be required in order to address either the contribution of mosaic or noncoding variation in these diseases.

uORF-introducing variants in the 5'UTR of the NIPBL gene as a cause of Cornelia de Lange syndrome.

Cornelia de Lange syndrome (CdLS) is a clinically-recognizable rare developmental disorder. About 70% of patients carry a missense or loss-of-function pathogenic variant in the NIPBL gene. We hypothesized that some variants in the 5'-untranslated region (UTR) of NIPBL may create an upstream open reading frame (uORF), putatively leading to a loss of function. We searched for NIPBL 5'-UTR variants potentially introducing uORF by (i) reannotating NGS data of 102 unsolved CdLS patients and (ii) literature and variant databases search. We set up a green fluorescent protein (GFP) reporter assay and studied NIPBL expression in a lymphoblastoid cell line (LCL). We identified two variants introducing a novel ATG codon sequence in the 5'-UTR of NIPBL, both predicted to introduce uORF: a novel c.-457_-456delinsAT de novo mutation in a 15-year-old male with classic CdLS, and a c.-94C>T variant in a published family. Our reporter assay showed a significant decrease of GFP levels in both mutant contexts, with similar levels of messenger RNA (mRNA) as compared to wt constructs. Assessment of LCL of one patient showed consistent results with decreased NIPBL protein and unchanged mRNA levels. 5'-UTR uORF-introducing NIPBL variants may represent a rare source of pathogenic variants in unsolved CdLS patients.

Disentangling molecular and clinical stratification patterns in beta-galactosidase deficiency.

INTRODUCTION: This study aims to define the phenotypic and molecular spectrum of the two clinical forms of ß-galactosidase (ß-GAL) deficiency, GM1-gangliosidosis and mucopolysaccharidosis IVB (Morquio disease type B, MPSIVB). METHODS: Clinical and genetic data of 52 probands, 47 patients with GM1-gangliosidosis and 5 patients with MPSIVB were analysed. RESULTS: The clinical presentations in patients with GM1-gangliosidosis are consistent with a phenotypic continuum ranging from a severe antenatal form with hydrops fetalis to an adult form with an extrapyramidal syndrome. Molecular studies evidenced 47 variants located throughout the sequence of the GLB1 gene, in all exons except 7, 11 and 12. Eighteen novel variants (15 substitutions and 3 deletions) were identified. Several variants were linked specifically to early-onset GM1-gangliosidosis, late-onset GM1-gangliosidosis or MPSIVB phenotypes. This integrative molecular and clinical stratification suggests a variant-driven patient assignment to a given clinical and severity group. CONCLUSION: This study reports one of the largest series of b-GAL deficiency with an integrative patient stratification combining molecular and clinical features. This work contributes to expand the community knowledge regarding the molecular and clinical landscapes of b-GAL deficiency for a better patient management.

Neuropathological hallmarks of fetal hydrocephalus linked to CCDC88C pathogenic variants.

The prevalence of congenital hydrocephalus has been estimated at 1.1 per 1000 infants when including cases diagnosed before 1 year of age after exclusion of neural tube defects. Classification criteria are based either on CSF dynamics, pathophysiological mechanisms or associated lesions. Whereas inherited syndromic hydrocephalus has been associated with more than 100 disease-causing genes, only four genes are currently known to be linked to congenital hydrocephalus either isolated or as a major clinical feature: L1CAM, AP1S2, MPDZ and CCDC88C. In the past 10 years, pathogenic variants in CCDC88C have been documented but the neuropathology remains virtually unknown. We report the neuropathology of two foetuses from one family harbouring two novel compound heterozygous pathogenic variants in the CCDC88C gene: a maternally inherited indel in exon 22, c.3807_3809delinsACCT;p.(Gly1270Profs*53) and a paternally inherited deletion of exon 23, c.3967-?_c.4112-?;p.(Leu1323Argfs*10). Medical termination of pregnancy was performed at 18 and 23 weeks of gestation for severe bilateral ventriculomegaly. In both fetuses, brain lesions consisted of multifocal atresia-forking along the aqueduct of Sylvius and the central canal of the medulla, periventricular neuronal heterotopias and choroid plexus hydrops. The second fetus also presented lumbar myelomeningocele, left diaphragmatic hernia and bilateral renal agenesis. CCDC88C encodes the protein DAPLE which contributes to ependymal cell planar polarity by inhibiting the non-canonical Wnt signaling pathway and interacts with MPDZ and PARD3. Interestingly, heterozygous variants in PARD3 result in neural tube defects by defective tight junction formation and polarization process of the neuroepithelium. Besides, during organ formation Wnt signalling is a prerequisite for planar cell polarity pathway activation, and mutations in planar cell polarity genes lead to heart, lung and kidney malformations. Hence, candidate variants in CCDC88C should be carefully considered whether brain lesions are isolated or associated with malformations suspected to result from disorders of planar cell polarity.

NGLY1 Deficiency: A Rare Newly Described Condition with a Typical Presentation.

NGLY1 deficiency is the first recognized autosomal recessive disorder of N-linked deglycosylation (NGLY1-CDDG). This severe multisystemic disease is still poorly known and, to date, most cases have been diagnosed through whole exome or genome sequencing. The aim of this study is to provide the clinical, biochemical and molecular description of the first NGLY1-CDDG patient from France along with a literature review. The index case presented with developmental delay, acquired microcephaly, hypotonia, alacrimia, feeding difficulty, and dysmorphic features. Given the complex clinical picture and the multisystemic involvement, a trio-based exome sequencing was conducted and urine oligosaccharides were assessed using mass spectrometry. The exome sequencing revealed a novel variant in the NGLY1 gene in a homozygous state. NGLY1 deficiency was confirmed by the identification of the Neu5Ac1Hex1GlcNAc1-Asn oligosaccharide in the urine of the patient. Literature review revealed the association of some key clinical and biological features such as global developmental delay-hypertransaminasemia, movement disorders, feeding difficulties and alacrima/hypolacrima.

Hypersociability associated with developmental delay, macrocephaly and facial dysmorphism points to CHD3 mutations.

CHD3-related syndrome, also known as Snijders Blok-Campeau syndrome, is a rare developmental disorder described in 2018, caused by de novo pathogenic variants in the CHD3 gene. This syndrome is characterized by global developmental delay, speech delay, intellectual disability, hypotonia and behavioral disorders including autism spectrum disorder (ASD). Typical dysmorphic features include macrocephaly, hypertelorism, enophthalmia, sparse eyebrows, bulging forehead, midface hypoplasia, prominent nose and pointed chin. To our knowledge, there have been no other clinical descriptions of patients since the initial publication. We report the clinical description of a 21-year-old patient harboring a pathogenic de novo variant in CHD3. We reviewed the clinical features of the 35 previously reported patients. Main features were severe intellectual disability, dysmorphic facies, macrocephaly, cryptorchidism, pectus carinatum, severe ophthalmologic abnormalities and behavioral disorders including ASD, and a frank happy demeanor. Hypersociability, which was a noticeable clinical feature in our case, despite ASD, is an uncommon behavioral feature in syndromic intellectual disabilities. Our report supports hypersociability as a suggestive feature of CHD3-related syndrome along with developmental delay, macrocephaly and a dysmorphic facies.

X-linked partial corpus callosum agenesis with mild intellectual disability: identification of a novel L1CAM pathogenic variant.

Pathogenic variants in L1CAM, the gene encoding the L1 cell adhesion molecule, are responsible for a wide clinical spectrum including X-linked hydrocephalus with stenosis of the Sylvius aqueduct, MASA syndrome (mental retardation, aphasia, shuffling gait, adducted thumbs), and a form of spastic paraplegia (SPG1). A moderate phenotype with mild intellectual disability (ID) and X-linked partial corpus callosum agenesis (CCA) has only been related to L1CAM in one family. We report here a second family, including 5 patients with mild to moderate ID and partial CCA without signs usually associated with L1CAM pathogenic variations (such as hydrocephalus, pyramidal syndrome, thumb adductus, aphasia). We identified a previously unreported c.3226A > C transversion leading to a p.Thr1076Pro amino acid substitution in the fifth fibronectin type III domain (FnIII) of the protein which co-segregates with the phenotype within the family. We performed in vitro assays to assess the pathogenic status of this variation. First, the expression of the novel p.Thr1076Pro mutant in COS7 cells resulted in endoplasmic reticulum (ER) retention and reduced L1CAM cell surface expression, which is expected to affect both L1CAM-mediated cell-cell adhesion and neurite growth. Second, immunoblotting techniques showed that the immature form of the L1CAM protein was increased, indicating that this variation led to a lack of maturation of the protein. ID associated with CCA is not a common clinical presentation of L1CAM pathogenic variants. Genome-wide analyses will identify such variations and it is important to acknowledge this atypical phenotype.

Detection of copy-number variations from NGS data using read depth information: a diagnostic performance evaluation.

The detection of copy-number variations (CNVs) from NGS data is underexploited as chip-based or targeted techniques are still commonly used. We assessed the performances of a workflow centered on CANOES, a bioinformatics tool based on read depth information. We applied our workflow to gene panel (GP) and whole-exome sequencing (WES) data, and compared CNV calls to quantitative multiplex PCR of short fluorescent fragments (QMSPF) or array comparative genomic hybridization (aCGH) results. From GP data of 3776 samples, we reached an overall positive predictive value (PPV) of 87.8%. This dataset included a complete comprehensive QMPSF comparison of four genes (60 exons) on which we obtained 100% sensitivity and specificity. From WES data, we first compared 137 samples with aCGH and filtered comparable events (exonic CNVs encompassing enough aCGH probes) and obtained an 87.25% sensitivity. The overall PPV was 86.4% following the targeted confirmation of candidate CNVs from 1056 additional WES. In addition, our CANOES-centered workflow on WES data allowed the detection of CNVs with a resolution of single exons, allowing the detection of CNVs that were missed by aCGH. Overall, switching to an NGS-only approach should be cost-effective as it allows a reduction in overall costs together with likely stable diagnostic yields. Our bioinformatics pipeline is available at: https://gitlab.bioinfo-diag.fr/nc4gpm/canoes-centered-workflow .